RECURRENCE OF PRODUCT OF LINEAR RECURSIVE FUNCTIONS 

CHENG LIEN LANG AND MONG LUNG LANG 

Abstract. We study the recurrence of the product of n functions, each of which satisfies 
("^ , the recurrence relation x{m) — Aix{m — 1) + A2x{m — 2) + • • • + Asx{m — s). 

(N 

j>^! 1. Introduction 

Let {x{m) : tti G Z} be a function defined by the following linear recurrence relation 



C^ 



^O I x{m) = Aix{m — 1) + A2x{m — 2) + • • • + Asx{m — s), (1-1) 






where Ai,A2,--- ,As are constants and Ai ^ for some i. Denoted by {u{m) : m €z 7j} 
the function that satisfies (1.1) and that u{l) = u{2) = • • • = u{s — 1) = 0, u{s) = 1. The 
main purpose of this article is to determine a recurrence relation for X{m), where X{m) is a 

1 -Q , product of n functions, each of which satisfies (1.1). The closed form (in terms of matrix) of 

d ' such a recurrence relation for X{m) can be found in Section 2. This fact answers a question 

(and its generalisation) raised by Cooper and Kennedy [CK] which was partially solved by 
Stinchcombe [S] (Partially solved is the sense that the characteristic polynomial of his Xm 

f^ ' must admit distinct roots). Two amazing recurrences can be found in Section 5 of [S]. The 

^ . main result of our study can be found in Proposition 2.3. Our method is elementary and uses 

^'^ ' linear algebra only. 

oo ; B J- 

^O ! Section 3 studies the recurrence relation of the product of n recursive functions, each of 

f^ ' which satisfies the recurrence Wm = pWm-i — qWm-2- Note that such recurrence has been 

'nJ" . studied by Jarden [J], Brenann [B] and Cooper and Kennedy [CK]. The main results can 

be found in Proposition 3.2 and 3.2. Our proof is direct and is slightly different from [CK]. 
Product of two recursive sequences each of which satisfies Wm = pWm-i — qWm-2 is the most 
interesting case. The recurrence for such functions can be used to verify various identities 
concerning Fibonacci and Lucas numbers. See [LLl] and [LL2] for more detail. 

A rather interesting connection (from our point of view) between the identity L^ — 5F^ = 
^_' 4(— 1)" and the Galois field of the characteristic polynomial of the recurrence of X{m) = F^ 

can be found in Section 4. 

2. The matrix form of the recurrence of X{m) 

The main purpose of this section is to give a recurrence relation for X{m) (a product of n 
functions Xi{m), each of which satisfies (l-l))- We shall start with an example which best 
reveals our idea and strategy. 

Example 2.1. Suppose that x{m) satisfies the recurrence x{m + 2) = px{m + 1) — qx{m). 
Let X{m) = x(m)'^. Then X{m) satisfies the following recurrence relation. 

X{m + 3) = (/ - q)X{m + 2) + {q"^ - p^q)X{m + 1) + q^X{m). (2.1) 

Proof. Denoted by u{m) the recursive function that satisfies u{m + 2) = pu{m + 1) — qu(m) 
and ti(0) = 0, n(l) = 1. One sees easily that 

1 



O 

m 



X 



x{m + r + 1) = u{m + l)x{r + 1) — qu{m)x{r). (2-2) 

Take the square of the left and right hand side of (2.2) and shift all the squares to the left of 
the equation, one has 

x(r7T, + r + 1) — u{m + 1) x[r + 1) — q u{m) x{r) =2qu{m + l)u{m)x{r)x{r + 1). (2.3) 

We note that on the right hand side of (2.3), the function A{r) = 2qx{r)x{r + l) is independent 
of ra. The cases for m = 1 and 2 are given as follows. 

x(r + 2f - u{2fx{r + if - q^u{lfx{rf = 2qu{2)u{l)x{r)x{r + 1), (2.4) 

x{r + 3)2 - ui3fxir + 1)^ - q'^u{2fxirf = 2gu(3)u(2)x(r)x(r + 1). (2.5) 

We now put these two equations into the following 2x2 matrix. 



Z 



x{r + 2)2 - u{2fx{r + 1)^ - q'^u{lfx{r)'^ u{2)u{l) 
x(r + 3)2 - u{3fx{r + if - q^u{2fx{rf u(3)ii(2) 



(2.6) 



Denoted by Cj the columns of Z. Applying (2.4)and (2.5), the first column is a multiple of the 
second column. To be more accurate, Ci — 2qx{r + l)x{r)C2 = Ci — A{r)C2 = 0. Hence the 
determinant of the above matrix Z is zero. Note that this holds for all r and that u{l) = 1, 
u{2) = p and u(3) = p^ — q- This completes the study of our example. 

We shall now study the general case by a simple generalisation of the techniques we presented 
in Example 2.1. 

Lemma 2.2. Let x{m) he a function satisfies (1,1). Then 

x{rn + r + 1) = ai{m)x{r + 1) + a2{rn)x{r) + • • • + as{rn)x{r + 2 — s), (2-7) 

where ai{m) = u{rn + l),ai{m) = Aiu{'m) + Ai^iu{rn — 1) + • • • + Asu{rn — s + i) for i > 2. 
In particular, ai{m) satisfies (1.1) for every i. 

Proof. Apply mathematical induction. D 

Let xi{m),X2{m), ■ ■ ■ , Xnim) be functions, each of which satisfies (1.1). It follows from Lemma 
2.2 that each one of them satisfies (2.7). Let X{m) = xi{m)x2{m) ■ ■ •x„(?n). Applying (2.7), 
one has the following. 

s s 

X(m + r + 1) - ^ ai{mTX{r + 2 - z) = ^ A(ei, ea, • • • , e,) [] ai{mY\ (2.8) 

1=1 A i 

where A is the set of (ei, 62, • • • , es)'s such that < Cj for alH, 61 + 62 + - • - + 6^ = n and 6^/0 
for at least two i's as we have shifted all the X{r + 2 — i)'s to the left of the equation. A has 
k members (see Remark 2.4). ^(61, 62, • • • , Cg) is a function in Xj(r + j)'s, where 1 < i < n, 
2 — s < J < 1 (take X{rn) = x{rnf in Example 2.1 for example, A{r) = 2qx{r + l)x(r)). It is 
important to note that A{ei, 62, • • • , 6^) is independent of m. For instance, in (2.8), regardless 
the value of m, the coefficient of as-i{m)as{m)^~^ is given by x(r + 2 — (s — l))x(r + 2 — s)""-*^. 
Identity (2.8) for m = 1, 2, 3, • • • take the following forms 

2 



i^l A « 

s s 

X{2 + r + l)-J2 a^{'^TX{r + 2 - z) = ^ A(ei, 62, • • • , e,) [] a,(2)^* (2.82) 

i^l A i 

s s 

X(3 + r + 1) - ^ a.(3)"X(r + 2 - 2) = ^ A(ei, 62, • • • , e,) [| a,(3)^- (2.83) 



Set 



i=l 



c(m) = X{m + r + 1) - ^ ai{mYX{r + 2-i). 

i=l 



(2.9) 



Associate to (2.8m) a row vector of length A; + 1 of the foUowing form (the terms Y\1 ai{mY^ 
are ordered according to the lexicographical order of A) 

Vrn = (c(m),ai(m)""^a2(m), ,as^2{rn)as{mY^^ ,as-i{m)as{m)'^^'^). (2.10) 

Similar to how we form matrix (2.6) from (2.4) and (2.5), we now put these k + 1 vectors 
^1)^1+1) • • • j^'i+fc (in this order) into the following (A; + 1) x (A; + 1) square matrix which we 
call it the recurrence matrix associated to Xim). 



Z 



c(l) 
c(l + l) 

c{l + k) 



a,_2(l)a,(l)"-i a,_i(l)a,(l)-i 

a,_2(l + l)a,(l + 1)^-1 a,_i(l + l)a,(l + l^-^ 



as-2{l + k)as{l + ky-^ as-i{l + A;)a^(l + k) 



n-l 



(2.111 



Proposition 2.3. Let X{m) be a product of n functions, each of which satisfies (1.1). Then 
X{m) admits a recurrence relation. Let Y{m) he another function which is also a product of 
n functions, each of which satisfies (1, 1). Then X{m) and Y{m) satisfy the same recurrence 
relation. 

Proof. Denoted by Ci the first column of Z and denoted by Ceie2---es the remaining columns. 
Applying identities (2.8i)-(2.8fc_|-i) and (2.9), the first column of the above matrix Z is a linear 
combination of the remaining columns. To be more accurate, 

Cl - ^^^(61,62, 



A 



I CsjL^eie2---es — U. 



2.12 



As a consequence, the determinant of Z is zero. Hence the set of row vectors {-Ri, i?2) • " " j-^fc+i} 
of Z is a linearly dependent set. Hence there exists Tj {0 < j < k), not all zero, such that^ tj-Rj 
takes the form (X^.-^g '''j'^i^ + i)' 0) 0) " " " > 0)- Since the fact that the first column is a linear 
combination of the remaining columns remains valid under row operation, we conclude that 
^rjc(l + j) = 0. We note that in (2.7)-(2.9), r is independent of m. This implies that 
= E^ic(l+j) = E'ri[^(l + ^ + l)-Ei=i«i(l)"-^(^ + 2-?)] for ah r. As a consequence, 
this gives a recurrence relation for X{m). 

Since X(m) and Y{m) admit the same recurrence matrix Z, they satisfy the same recurrence 
relation. D 



Remark 2.4. Apply the Multinomial Theorem for (yi + 2/2 + • • • Us)^, k 



n + s - 



Remark 2.5. (i) Denoted by Ci, C2, • • • , Ck+i the column vectors of Z. Let V = [C2, C3, • • • , 
Ck+i]- In the case the matrix V = [C2,C3,--- ,Ck+i] is of rank k, the cofactor expansion 
along the first column of Z gives a recurrence relation of X{m). In the case the rank of V is 
r and r < /c — 1, one can still apply cofactor expansion to an (r + 1) x (r + 1) submatrix of Z 
to get a recurrence relation of X{m). 

(ii) In [S], Stinchcombe [S] studies the characteristic polynomial of (1.1) and draws the conclu- 
sion about the existence of a recurrence relation for X{m) = x(m)". Note that his method is 
different from ours where his calculation must be carried through over the complex field C and 
the roots of the characteristic polynomial of x{m) must be distinct. Two amazing recurrences 
can be found in Section 5 of [S]. 

2.1. Discussion, (i) Let X{m) be given as in Proposition 2.3. Suppose that the recurrence of 
X{m) involves t terms. Since {X(l + /c), • • • , X(t+/c)} satisfies the same recurrence relation for 
arbitrary fe, determinant of the following matrix is zero, which can be viewed as a recurrence 
relation for X{m) as well. This generalises a result stated in [W]. 



X{k + l) X{k + 2) ■■■ X{k + t) 

X{k + t + l) X{k + t + 2) ■■■ X{k + 2t) 

X{k + f -t + l) X{k + f -t + 2) ■■■ X{k + t^) 



(2.13) 



(ii) As one recurrence function may satisfy more than one recurrence relation, one would like 
to know whether Proposition 2.3 offers the best recurrence relation for X{m). The answer is 
No. Take x{m) = 2x{m — 1) — x{m — 2) {x{Q) = x{l) = 1) for instance, applying Proposition 
2.3, the recurrence for X{m) = x(m)^ is X{m) = 3X(m — 1) — 2>X{m — 2) + X{m — 2). 
This is obviously not the best recurrence for X{m) as X{m) satisfies the recurrence relation 
X{m) = 2X{m — 1) — X{m — 2) as well. Is this because the characteristic polynomial of 
x{m) = 2x{m — 1) — x{m — 2) has repeated root ? Again, the answer is No. Take y{m) = 
2y{m — 1) — y{rn — 2) (y(0) = 0, |/(1) = 1) for instance, Y{rn) = y{rn)'^ satisfies the relation 
Y{m) = 3y (m - 1) - ?,Y{m - 2) + Y{m - 2) but Y{m) / 2Y{m) - Y{m - 1). Hence the initial 
values also play some roles. 

3. Recurrence Relation for powers of Wm = pWm~-i - qWm-2 

Let a,b,p,q £ C, q ^ 0. Following the notations of [Ho], we define the generalised Fibonacci 
sequence {Wm} = {Wmia, b;p, q)} by Wq = a, Wi = b, 

Wm = pWm-1 - qWm-2- (3.1) 

Obviously the definition can be extended to negative subscripts ; that is, for n = 1, 2, 3, • • • , 
define 

W-m = {pW-m+1 - W-m+2)/q. (3.2) 

In the case a = 0, b = 1, we shall denote the sequence {Wm{0, 1 ; p, q)} by {um}- Equivalently, 

Um = Wm{0,l;p,q). (3.3) 



Note that u. 
that 



-q "^Um (this simple fact can be proved by mathematical induction) and 



W, 



m+r+l 



^m+lWr+l 



qUmWr- 



(3.4) 



The purpose of this section is to give a recurrence relation of X{m), a function which is a 
product of n functions, each of which satisfies the recurrence (3.1). Applying our results in 
Proposition 2.3, the recurrence matrix Z (see (2.11)) is an n x n matrix of the following form 
and detZ = gives us a recurrence for X{m). A detailed study of the recurrence detZ = 
can be found in sections 3.1 and 3.2. 



Z 



X{r + 2) - Xir + l)^^ - i-q)''X{r)u1 
X{r + 3) - X{r + 1)^' - {-q)''X{r)u^ 



"2 ^«i 



ri-2^,2 



X{n + r + l)- X{r + l)<+i - (-q)"X(r)< u. 



"-1,, ,,"-2 2 

n+l"n "n+l"n 



13-^2 



^n+l"n 



(3.5) 



3.1 The nontrivial case. Throughout this subsection, we shall assume that U1U2 ■ ■ ■ "Un+i 7^ 0. 
Denoted by Cj, the i-th column of the above matrix. The determinant of Z can be determined 
as follows. Let R = (n^,u^,--- ,<+i)*, S = (ny,^^,--- ,<)*, T = (A:(r + 2),X(r + 3), • • • , 
X{n + r + 1)Y (t for transpose). Set A = (R, 02,03, •••c„), B = (S*, C2, C3, • • • c„), C = 
(T, C2, C3, • • • On)- Since the determinant function is a linear function, one has 

det Z = detC - X{r + 1) • det^ - {-qfX{r) • detB = 0. (3.6) 

Applying {A3) of Appendix A, Appendix B and (C2),(C5) of Appendix C, one has the 
following. 

Lemma 3.1. Let n G N 6e fixed and let A, B and C be n x n matrices given as above {see 
(3.5) and (3.6)). Suppose that U1U2 ■ ■ ■ ^in+i 7^ 0. Then 



detC = ^Ci, detB = {-l)"-~^a{n-l) det Ai, detA = a{n)detAi, 



(3.7) 



j=i 



where Ai is a nonsingular n x n matrix, cr{n) = Un+iUn ■ ■ ■ U3U2 and the quantity Ci is given 
by the following. 



Ci 



X{r + i + l)a{n)a{n - l)detAidetQ\^^ 



(3.8) 



Ui+iUi{ui-iU2-i • • • U^i){uiU2 ■ ■ ■ Un-i)detQ\^ ' 

where 1 < i < n, u_m = ~Q~"^u„i {q is the fixed constant given in (3.1)) and Q„i is an m x m 
matrix of determinant g™(™-i)/2^ 

Proof. See Appendix A, B and C. D 

Note that detyli appears in Ci, det-B and det^ and that detAi is nonzero. Applying (3.7), 
(3.8) and u-m = —q~"^Um, one may now simplify (3.6) and conclude that 

n+l 

J2^_lyqi(^-l)/2(^^ + l\i)uX{n + r + l-i)=0, (3.9) 

i=0 

where {m\k)u = 1 if /c = and {m\k)u = UmUm-i ■ ■ ■ Um-k-i / UkUk-i ■ ■ ■ ui \i \ < k < m. 
Note that {m\k)u is known as the generalised binomial coefficient. In summary, the following 
proposition holds. 



Proposition 3.2. Let n &N be fixed and let X{m) be a product of n functions, each of which 
satisfies (3.1). Suppose that uiU2 ■ ■ ■ lin+i ^ 0. Then X{m) satisfies the following recurrence. 

n+l 

^(_l)igi(i-i)/2(^ + l\i)uX{m - i) = 0. (3.10) 



Proof. Alternative Proof of Proposition 3.2. One sees easily that the X{m) in Proposition 
3.2 and ■uJJi admit the same recurrence matrix (see (2.11) for the definition and (3.5) for the 
actual matrix). Hence X{m) and uj^ satisfy the same recurrence relation (see Proposition 2.3). 
One may now prove the proposition by applying Jarden's Theorem which gives the recurrence 
for u^ (see Theorem 4.1 of [CK]). D 

3.2. The trivial case. Throughout this section, uiU2 ■ ■ ■ Un+i = 0. This implies that u^ = 0, 
for some k, where 1 < A; < n+l. Identity (3.4) now becomes Wk+r+i = Uk+iWr+i- Since X{m) 
is a product of n function, each of which satisfies Wk+r+i = Uk+iWr+i, one has X{k + r+l) = 
u^_^-^^X (r + 1) . Since this holds for every r G Z and k is fixed, the following proposition is clear. 

Proposition 3.3. Let n £N be fixed and let X{m) be a product of n functions, each of which 
satisfies (3.1). Suppose that Uk = 0, for some k, where 1 < k < n + 1. Then X{m) satisfies 
the following recurrence. 

X{k + r+l)=ul^^Xir + l). (3.11) 

Example 3.4. Let uq = l,ui = 1 and n„ = 2un-i — 4:Un-2- Then U2 = 2,U3 = 0, n4 = —8. 
By Proposition 3.3, u^ satisfies the following recurrence 

tx^+4 + 8^^r+i = 0. (3.12) 

Note that u^ also satisfies the recurrence relation n^+4 — 8n^_,_3 — 512ti^_^;^ + 4096n^ = 0. 

3.3. Application. The following corollary of Proposition 3.2 has been applied by Lang 
and Lang ([LLl] and [LL2]) to prove various identities concerning the generalised Fibonacci 
sequence. 

Corollary 3.5. A{m) = W2m., B{m) = WmWm+r and C{m) = q™- satisfy the following 
recurrence relation 

X{m + 3) = (p2 - q)X{m + 2) + (q^ - p^q)X{m + 1) + q^X{m). (3.13) 

The following corollary of Proposition 3.2 shows that the recurrence (3.10) can be used 
to describe the characteristic polynomial of Qv (see (^4) of Appendix A) which reveals the 
fact that the recurrence relation carries a lot of information about the generalised Fibonacci 
numbers. Proof of Corollary 3.6 can be found in Appendix D. 

Corollary 3.6. Let Q^ be given as in (Ai) of Appendix A. Then the characteristic polynomial 
of Q-o is X]i'=o(~^)*9* ("^N)"-^*- -^ote that f{x) is the polynomial that characterises the 

recurrence (3.10). 



4. The Galois Group of the Recurrence Relation 
Let p and q be given as in (3.1) and let 



rip,q,x) = ^i-iyq'^'~'^/Hn\i)ux\ (4.1) 



j=0 

We call the polynomial in (4.1) (see Corollary 3.6 also) the Galois polynomial of the recurrence 
relation (3.10). Similarly we call the Galois group Gnip,q) of (j)nip,Q,x) over Q the Galois 
group of the recurrence relation (3.10). 

Proposition 4.1. G„(l, —1) — Z2, a cyclic group of order 2, for all n > 2. 

Proof. Let a and h be roots of (/)2(x) = x^ — x — 1 = 0. One can show by induction that (see 

[B]) 

0,(1, -1, x) = {abr-\x - a")(x - b^)<Pn-2{l, -1, x/ab) (4.2) 

for all n > 2 {(l)o{x) is defined to be 1). Hence (/>n(l, —1, x) splits completely in the field Q(v5) 
andG„(l,-l) ^^2. D 

4.1. Discussion. Let -F„ and L„ be the n-th Fibonacci and Lucas numbers respectively. It is 
a well known fact that 

Ll-5F^=4{-ir. (4.3) 

Identity (4.3), from our point of view, is a must rather than an interesting connection between 
Fibonacci and Lucas numbers as the following suggested. Recall another well known fact about 
the factorisation of (pni^i — l,^;)- 

Mh -1, X) = (-l)"(x2 - LnX + (-1)")0„_2(1, -1, -X). (4.4) 



Since 0„(1, -l,x) splits in Q(\/5) and x^ - LnX + (-1)" = splits in Q{^/Ll - Ai-l)"^), one 
must have 

Ll-A{-ir = 5Al (4.5) 

for some An G N. This tells us that the difference between 4(— 1)" and the square of the 
Lucas number L„ must be five times a square A^. As for why An must be -F„, we note that 
both Fn and L^ satisfy the recurrence relation (3.13) (to be more accurate, the recurrence 
X{n + 3) = 2X{n + 2) + 2X{n + 2) - X{n)) and that L^ - 4(-l)" and F^ have the same 
initial values. In general, one has (Lemma 3.3 of Cooper and Kennedy [CK]), 

n 

MP,Q,x) = ll{x-a^f3^-^), (4.6) 

i=o 

where a and /? are roots oi x'^ — px + q = 0. Hence G{p, q) = Z,2 for all n > 2. This fact will 
give an identity similar to (4.5). 



5. Appendix A 

Throughout the appendix uiU2 ■ ■ ■ Un+i 7^ 0. We shall give full detail of how the determinant 
of the matrix A is evaluated {A is given as follows). Applying our technique given in this 
appendix, the determinants of B and C can be calculated similarly (see Appendix B and C). 



A 






.n-l„ 



'U2 U^ 






Ur,U 



n— 1 



2"-! 



.n-l„ 



n—2 2 



,n— 1 



,n-l 



(Al) 



Note that Uj+i is a common factor of the entries of the i-th row. Hence det A can be written 
as the product Un+iUn • • • U2 det Ai, where 



A^ 






Uo 



2 "^2 ""l 

n—1 „,n— 2„ 



ur^n? 



"2 



Uo U2 Un 



n-1 
L "ra+l 



U. 



U„,i u. 



n+l"'n "'n+l""n 



U„,-i Ur, 



n-1 
1 
,n— 1 



U 

Un 



u: 



n—l 



iA2) 



For our convenience, we shall define cr{n), T{n) as follows which will be used in the following 
discussion. 

n+l 

a{n) = Ylui, r(n) = ^"("-1)/^ det A = a{n) det Ai. {A3) 

i=2 

The rest of this section is devoted to the determination of the determinant of Ai. To save 
space, we denote the binomial coefficient ( ^ ) by (n|A;). Consider the following matrix. 



Qr. 



(n-l|l)p«-2(-g) 



(n-2|lK-^(-(7) 
(n-2|2)p"-4(_g)2 



(n - l|n - 2)p(-g)"-2 („ _ 2|n - 2)(-g)"-2 
(n-l|n-l)(-g)"-i 



P 1 

-q 










(A4) 



Note that Qn is an ra x n matrix where the i-th column of Qn gives the coefficients of the 
binomial expansion of [px — qy)"'~^ and that det Qn = T{n) = gr"-("--i)/2^ Since Ur+i = 
pUr — qur-i, uq = 0, Ul = 1 (scc (1.3)), the multiplication of Ai by Qn shifts the indices. 
That is. 



AiQ-' = 



.n-1 



,n-2„ 



Ul 



^rt ^n ^-r 



-'n-l 



AlQn = 



^n+2 



,n-2„ 



"3 



ri-2 
''^n+2^n+l 



-^n+1 



(A5) 



Since uq = 0,ui = 1, the first row of AiQ^^ is (l"^-*^, 0, 0, • • • ,0). Applying tlie cofactor 
expansion to tlie first row of (A5), the determinant of AiQ~^ is the determinant of the following 
matrix. 



Z 






1 

'U2 



U„ 



^n-1 



«2 ^uf 



Un 



u 



w: 



-3, ,2 



*n-l 






n— 1 



(^6) 



Compare the matrices in (^1) (an n x n matrix) and (^6) (an (n — 1) x (n — 1) matrix), 
we have established a recursive process which enables us to calculate the determinant of the 
matrix A as well as Ai. 



det Ai = T{n)T{n -!)••• T{2)a{n - 2)a{n - 3) • • • a{l) / 0. 



(^7) 



Since the multiplication of Ai by Qn shifts the indices (see (^5)), the following is clear. Note 
that this matrix plays an important in the study of matrix C (see Appendix C). 



AiQ; 



ra-l 


n-2 


ra-3 2 


«2-i 


U2-i ^1-i 


^2-i ^I 


n-1 


n-2 


n-3 2 


%-i 


%_j U2-i 


^3-i ^2 



,n-l 



n— 1 



L «n+l- 



U, 



^ra+1- 



-=^n2 



^i; 



n— 3 



.•^i„ 



,n-l 

n-1 
2-i 






u, 



n-1 



U, 



n-1 



(A8) 



Note that the i-th row of AiQ^"^ takes the form (1" , 0, 0, • • • ,0). The determinant of AiQ 
can be determined by the cofactor expansion by the i-th row (see (C5) of Appendix C). 

6. Appendix B 
Throughout the appendix uiU2 ■ ■ ■ "Un+i 7^ 0. The matrix B is given as follows. 



B 



u'{ 



.n-l„ 






U2 Ug 



-'n+l^n 






-2„,2 



n- 



UoU 



-1 1 



Uou: 



3 "-2 



n, 



n+l "n 



ti, 



n-1 



{Bl) 



After the removal of the common factor Ui from each entry of the i-th row, one sees that the 
resulting matrix is just the matrix Ai (see equation (^2)) when the first column is moved to 



the last. Hence deti? = (— 1) 



n-1 



a[n 



l)det ^1 (see A(3) for the definition of a{n — 1)). 



7. Appendix C 

Throughout the appendix uiU2 ■ ■ ■ Un+i / 0. The matrix C is given as fohows. Its deter- 
minant can be calculated by the cofactor expansion of the first column. The calculation is 
tedious but elementary. 



C 



W, 



r+2 



w. 



r+3 



^2 U^ U.^ 



.n-l„ 



™-2„,2 



-2„,2 



,n-2^,2 



"n+r+1 ^n+l'^^n ^n+l'^^n 



UoU 



n-1 



2"'l 



UnU. 



.n-1 



3 "2 



''^ra+l^n 



(CI) 



Denoted by Cj the i-th cofactor of the first column of C. The determinant of C is given by 
the following. 



det C7 = Ci + C2 + • • • + C„ 



{C2) 



where Ci = (— 1)*~'"^ W^j_,_]^det X and X is the (n — 1) x (n — 1) matrix that takes the following 
form. 



X ^ 



U2 



l-i+2 "i+1 



"n+l"n 






-^1+2 "i+1 



_1 tA^ 



U2<"' 



M,-U,- 



-'i+2"i+l 



-*n+l"n 



(C3) 



Note that Uj+iUj is a common factor of the entries of the j-th row when j < i — 1 and that 
Uj+2Uj+i is a common factor of the entries of the j'-th row when j > i. As a consequence, the 
determinant of X is the product of a{n)a{n — l)/ui-^-iUi and the determinant of Y, where 



Y 



^i+2 



n — 3^ 



ur^M.- 



71 — 3 ^ 



-'i+2 "i+1 



-1 "n+l"Ti 



n-2 
n-2 



YQ-U = 



-2 

hl-i 



:^, 



"0 



(C4) 



Similar to (A8), we consider the matrix YQ~^_^. It is an easy matter to write down the matrix 
as multiplication by Qn-i shifts the indices (see (^5) and (^8)). To calculate the determinant 
of y, we consider the cofactor expansion of the matrix AiQ~^ by the ith row (see (^8)), where 
the i-th row of AiQ~'^ takes the form (1"~^,0,0, • • • ,0). An easy observation of the actual 
forms of the matrices AiQ~'^ and YQ'Z'^ shows that 



det^iQ,/ = {-iy^^{ui^iU2-i---u^i){uiU2 



llfi 



,)detyQ-i,. 



(C5) 
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In summary, 



Ci 



W:;^i^^(j{n)(j{n - l)det Ajdet Qj^^ 

Ui+iUi{ui-iU2-i ■ ■ ■ U-i){uiU2 ■ ■ ■ Un-i)det Q\ ' 



(C6) 



Recall that u_m = —Q. "^Um (this simple fact can be proved by induction) and that detQn = 
^n(n-i)/2_ 'j^j^jg gives the determinant of C (see (C2)). 

8. Appendix D 

Suppose that uiU2 • • -u^ / 0. This appendix is devoted to the study of the characteristic 
polynomial of the matrix Q^ (see (^4)). For simplicity, we shall denote Qv by Q. Let E be 
the rational canonical form of f{x) = ^^"=0 ^i^^ — Z]r=o(~^)*9* (^l^)"-^*' Then 



E 






1 





•• 














1 


•• 

















1 •• 

















•• 





1 


ao 


Ol 


02 


03 •• 


• av-2 


av-1 



(Dl) 



Let Ai (set n = v) he the matrix given as in (^2). Applying Proposition 3.2, the function 
X{m) = u: 



??i— 1 "'m 



v—l—k 



satisfies the recurrence relation 



^i-i)y(^-^)/^(v\z)uul,ur'~' = o. 



{D2) 



1=0 



Applying {D2), the multiplication of Ai by E (to the left of ^i) shifts the indices. Since the 
multiplication of Ai by Q (to the right of Ai) shifts the indices as well (see (vl5)), one has 



EAi = AiQ, 



im 



which can be verified by direction calculation. Since the determinant of Ai is nonzero (see 
{A7)), E and Q are similar to each other. As a consequence, the characteristic polynomial of 
Q is /(x). 
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